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' We study the zero-temperature phase diagram of a square-lattice S = 1/2 Heisenberg anti- 

0^ , ferromagnet with two types of regularly distributed nearest-neighbor exchange constants, Ji > 

(antiferromagnetic) and —00 < J2 < 00, using spin-wave series based on appropriate mean-field 
Hamiltonian and exact-diagonalization data for small clusters. At a quasiclassical level, the model 
^-j, displays two critical points separating the Neel state from (i) a helicoidal magnetic phase for rela- 

tively small frustrating ferromagnetic couplings J2 < (J2/J1 < — 1/3 for classical spins), and (ii) 
a finite-gap quantum paramagnetic phase for large enough antiferromagnetic exchange constants 
J2 > 0. The quantum order-disorder transition (ii) is similar to the one recently studied in two- 
layer Heisenberg antiferromagnets and is a pure result of the zero-point spin fluctuations. On the 
other hand, the melting of the Neel state in the ferromagnetic region, J2 < 0, is a combined effect 
of the frustration and quantum spin fluctuations. The second-order spin-wave calculations of the 
ground-state energy and on-site magnetization are in accord with our exact diagonalization data 
in a range away from the quantum paramagnetic phase. In approaching the phase boundary, the 
■ theory fails due to the enhanced longitudinal spin fluctuations, as it has recently been argued by 

I ' Chubukov and Morr. 

m 

Q . PACS: 75.50.Ee, 75.10.Jm, 75.30.Kz, 75.10.-b 
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The zero-temperature properties of a number of two-dimensional Heisenberg spin models displaying order-disorder 
quantum phase transitions have recently attracted much attention. Most efforts were connected with the square- 
lattice Heisenberg antiferromagnets with frustrating second J 2 and third J3 nearest-neighbor bondstl as well as with 



a 

\ the two-layer quantum Heisenberg modeld 

The mentioned models present, concrete realizations of the quantum order-disorder transition considered by 
Chakravarty, Halperin, and Nelsoncl using a macroscopical approach based on the quantum non-linear a model in 2 + 1 
dimensions. Notice that in the J\ — J 2 — J3 models the destruction of the Neel state is a result of both the frustra- 
tion and zero-point spin fluctuations, whereas the two-layer antiferromagnet can be turned through an order-disorder 
transition by varying the antiferromagnetic exchange coupling between the planes. The latter model constitutes an 
ideal system where the zero-point spin fluctuations destroy the classical Neel state. n 

Less studied are the 2D Heisenberg antiferromagnets with random ferromagnetic exchange bondsa which are cur- 
rently considered as an appropriate model of the coppex-oxygen sheets of lightly doped insulating high-T c supercon- 
ductors. The idea comes from Aharony and co-workersQ who suggested that the holes introduced into the CuO layers 
are localized at individual oxygen sites. The coupling between the copper and oxygen spins results in an effective 
ferromagnetic coupling between the copper spins situated near an oxygen. 

In this paper we study a quantum Heisenberg antiferromagnet with two kinds of nearest-neighbor exchange cou- 
plings, regularly distributed on the square lattice, Fig.l. The relevance of the model is twofold. First, it combines the 
features of the J\ — J 2 — J3 and two-layer models in displaying the mentioned two mechanisms of destruction of the Neel 
phase and, on the other hand, the model can be considered as a prototype of some more realistic antiferromagnetic 
models containing ferromagnetic bonds. 
The Hamiltonian of the system reads 

H = J X H X + J 2 i/ 2 = Ji ^ SiSj + J 2 ^ S ; Sj , (1) 

<ij> (ij) 
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where the sums run over the nearest-neighbor J\ and J2 bonds, respectively. In what follows. we frequently use the 
notations a = J2/ Ji, Ji = 1- In the antiferromagnetic region J 2 > 0, Singh, Gelfand, and HuseQ have found for Eq.(jl]) 
a critical point a s ps 2.56 separating the classical Neel ordered and a finite-gap quantum paramagnetic phases, using 
series expansions around the dimer Hamiltonian H^,. They also argued that the critical exponents are consistent 
with those of the 3D classical Heisenberg model. In this study we analyze the zero-temperature phase diagram of 
the model in the range —00 < J2 < 00 both analytically, reducing the problem to an interacting Bose gas, and 
numerically, throughout exact diagonalization of small clusters. For ferromagnetic exchange bonds, at a = —1/3 for 
classical spins, we found another critical point, dividing the Neel phase from a helicoidal magnetic phase, Fig.l. The 
latter phase reduces to a three-sublattice Neel state in the limit J 2 = —00. The classical helicoidal phase is described 
by the equations 

S A (R) = u cos QR + v sin QR, 
S B (R + x) = ucos(QR + tt + 3$) +vsin(QR + 7r + 3$), (2) 

where the classical spins Sa and Sb belong to the sublattice A and B, respectively, u and v are perpendicular 
unit vectors in the spin space; R runs over the sites of the sublattice A. The pitch vector Q and the phase $ are 
determined through the relations 

cosQ x ( y ) = ^(-^--l\, Q y (x)=0, (3) 

t Qx(y) , 1 

The angle $ can be chosen as an order-parameter of the continuous classical transition. 

Below we concentrate on an analysis of the role of quantum spin fluctuations in the Neel phase of Eq.(Q). To this end, 
in the next section we follow a method overcoming the problems of the standard spin-wave theory (SWT) connected 
to logarithmic divergences at the classical transition points. Then we construct spin- wave series for the ground-state 
energy and on-site magnetization, and discuss the changes in the phase diagram implied by the interaction. In the 
last section we present exact diagonalization (ED) data for small clusters and compare them with the theory. 



II. SPIN- WAVE SERIES FOR THE GROUND-STATE ENERGY AND ON-SITE MAGNETIZATION 

We proceed with a short discussion of the results implied by the linear spin- wave theory (LSWT). LSWT predicts 
the existence of a region around the classical transition point a c = —1/3, where the magnetic states are completely 
melt by the zero-point spin fluctuations, Fig. 2. In this approximation, the spin-liquid phase exists for arbitrary S due 
to the logarithmic divergence at a — —1/3 in the integral 

5m = 1 — i — In 1 r , (4) 

N ^ e s k 2 \a — a c \ w 

defining the reduction of the on-site magnetization 5m = S— < S z >. Here N is the number of lattice sites, 
kq = (1 — a) /4. This behavior of 5m can be regarded as an aftereffect of the classical "soft lines" at the critical point 
in the low-energy magnon dispersion function 

e s k ~ Jcja^- a c )kl + c 2 k^, |k|«0,7r. (5) 

Due to the continuous character of the classical transition, one could expect the same behavior when the interaction 
is switched on. _A consideration valid in the large-S limit, and based on an appropriate quantum nonlinear a model on 
the critical lined, predicts the existence of a disordered spin-liquid phase in an exponentially small region around the 
critical line. There are some indications that the spin-liquid phase survives for arbitrary S up to the extreme quantum 
limit S = 1/23. In our case, however, the structure of the quasiparticle spectrum (i.e., the number of Goldstone 
modes) is changed, so that the problem requires a special studyu of the helicoidal phase, wich is beyond the scope of 
the present paper. 

LSWT also qualitatively predicts thalransition to a finite-gap quantum paramagnetic state in the antiferromagnetic, 
J2 > 0, region of the phase diagramEl However, the estimate a s sa 10.6 cpasiderably exceeds the series estimate 
a s pa 2.56. The situation is analogous to the one in two-layer antiferromagnetstil, i.e., SWT is not successful because of 
the enhanced longitudinal spin fluctuations. One of the purposes of the present study, using also exact diagonalization 
data, is to see to what extent SWT can be trust in approaching the quantum paramagnetic phase. Our model has 
the advantage of permitting exact diagonalization studies for large enough clusters. 
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A. The method 



As mentioned above, the standard spin- wave series are divergent at the classical transition point a c = —1/3. For 
example, the first-order correction to the on-site magnetization Am' 1 ' is expressed through the following divergent 
integrals 

(1 ) "o 1 {lAn. - «o cos fc^ (cos k x - f k ) - k sin 2 k x 
A - U -^^L i3 (6) 

k 



sk 



Kq v Kq(1 — Vk COS k x ) + ^k(^k — COS k x ) 



N ^ e s k 

k 



where = (cosfc x + cosfc y )/2, e s k = y(l — vE) — 2«o(l — ^k cos k x ). Since in the SWT formalism the on-site 
magnetization is a two-boson operator, it is clear that the divergences in the first-order correction Am' 1 ' come from 
two-boson interaction terms. A straightforward way to improve the situation is to use a normal-ordered operator 
form for the interaction. This approach leads to renormalization of the quasiparticle excitation spectrum and to some 
resumed SWT series starting now from an appropriate mean- field theory (MFT). Instead of 1/5, the new series are, 
however, governed by a formal small parameter A introduced in front of the quartic Dyson-Maleev normal-ordered 
boson interaction. Therefore, the quality of the series should be controlled through another method. In our study we 
present exact diagonalization data for small clusters. 

Let us briefly describe the procedure leading to the spin-wave series. SWT requires three subsequent transfor- 
mations. Through a standard Dyson-Maleev bosonization scheme, the original spin Hamiltonian transforms into 
equivalent boson Hamiltonian. In the usual SWT the two-boson terms scaling with 5 constitute the free-boson 
Hamiltonian Hq which fixes the 1/5 series. As it was mentioned above, in our case more convenient is to diagonalize 
the quadratic paxt-of the Hamiltonian by means of a generalized Bogoliubov transformation, also touching on the an- 
harmonic termsEjlij. The new quadratic Hamiltonian Ho can be obtained in the following way. First, we subsequently 
do the Fourier and quasiparticle Bogoliubov transformations, the latter being defined through the relations 

a k = Uka k + v k f3l, b k = u k /?k + fk4. l M k| 2 - \vk\ 2 = 1- (7) 

In a quasiparticle representation, the boson Hamiltonian takes the form 

H = E +H + Vdm, (8) 

where the Dyson-Maleev interaction Vdm is i n a normal-ordered operator form (see Appendix A). Then E is the 
mean-field ground-state energy and Hq is a quadratic free-particle Hamiltonian, containing also contributions from 
the operator ordering of the quartic interaction term. 

The diagonalization of Ho, which in terms of quasiparticle operators reads 



H ° = [#k(4«k + + ( fl ka k 4 + h.c.)] , (9) 

k 



requires i?k = 0. In principle, the condition Bk = leads to some integral equations for uit-,and «k studied in the 
general case by ObukhovE3. Recently, similar method has been applied to the J\ — Ji modeltj. 

Note that the so-defined mean-field theory can also be regarded as a result of the variational equations 

6 <m\*> j <w> = 

' \ „ ._ ^ iTf iTf v ' 



Su k < > ' Sv k < > 

(under the constraint |wk| 2 — \ v k\ 2 — 1), defining the generalized Bogoliubov transformation. The quality of the 
discussed perturbation scheme mainly depends on the quasiparticle interaction strength playing the role of a physical 
small parameter. Particularly, the latter is expected to increase near the critical points. The model Hamiltonian 
(|l|) represents a convenient system to check the theory due to the existing two critical points and a numerical-series 
estimated for the order-disorder transition point. 

B. Zeroth-order mean- field approximation 

The above program is straightforward for Eq.(|l|) and leads to a renormalization of the bare coupling (1 — J2/J1) — > 
(1 — J2/ J\)Z. The implicit equation for the factor Z = Z(S, a) reads 
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where 



a i 



1 . 1 1 - vv. cos k x 1 1 /ia\ 

2 + ' a 3 = --(l-«) + «a 1 + -^e k , (12) 

k k 

/ ^ 

£k = y (1 - ^) - 2k(1 - t'kcosfc^), K = (l-a) — . (13) 

The renormalization factor Z defines (up to an irrelevant phase factor) the Bogoliubov coefficients itk and = —v^/uk 



ll-K + Zk 
^ k -KCOSfc x 

3?(xk) = 7: — r, (14 

(1-K + Ek) 



(1-K + ek) 1 

Notice that one cannot choose the phase factor so as to make both u k and u k real functions. 
The quasiparticle excitation spectrum reads 



Eu = 4S (l + ^) e k , (15) 

where r is a numerical factor analogous to Oguchi's correction^, r = —2a 2 . 

The ground-state energy Eq and on-site magnetization mo in this approximation are 

|l = -(l-«Z)(S -*,)», (16) 

-^r^I? (17) 

k 

We will see below that the mean-field estimate for the ground-state energy E is in agreement with ED in a large 
region of the phase diagram. The on-site magnetization m vanishes at a c ~ —3.0, i.e., there is a shift of the classical 
transition point towards the helicoidal phase. This phenomenon is known as "Villain's order from disorder" fca. The 
transition to quantum paramagnetic state in MFT is located at a s « 5.7, which is approximately twice larger than 
the series estimate. Below we show that the second-order correction to this result improves the estimate. 



C. Second-order corrections to Eq and mo 

Now we proceed with a calculation of the second-order 0(V^ M ) corrections to E and tuq. It is easy to see that the 
first-order corrections both to Eq and mo identically vanish for the normal-ordered interaction Vdm- The second-order 
correction to Eq has the standard form (Appendix B) 

AEW /2\ 3 ^ a , (1234)T/( 8 ) (3412) , s 

= -( — > A (1 + 2-3-4) -, (18) 

2N \NJ f^ 4 V ' E 1 +E 2 +E- i + E i ' y ' 

where Ek are the quasiparticle energies Eq.@). The expression for the second-order correction to the on-site mag- 
netization consists of two parts (see Appendix B) 

Ato (2) = Ato^ 2) + Ato[ 2) , (19) 

where 
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a (2) 645,, r A(l + 2-3-4) 

x^ 7 )(1234)^ 8 >(3412) f i- + i- + i- + -M , (20) 



1-4 



W(1234) = y (7) (1234) 
+ F (8) (3412) 



5 )(3412) + 4^1 /(2) (3412) 



E2 E3 



y(6) (1234) + ?^iy(3) (1234) 



S4 



(22) 



In the last equations the sums run over the vectors ki , k2 , k3 , IC4 from the magnetic Brillouin zone; A(k) is the 
Kronecker function. The explicit expressions for the vertex functions = (k^^, kaki), i = 1, 9, are given in 
Appendix A. A straightforward task is to show that,the above 6D integrals defining AE^ and Am' 2 ' are convergent, 
using the vertex expansions near k = and k = mL 3 ]. 

The 6D integrals written above were calculated using a standard Gaussian method for a set of lattices with up to 
M = (20 x 20) 3 cells (each containing 4 points) and then the results were extrapolated to M — 00. Starting from the 
(4 x 4) 3 -cell lattice the results change monotonously and ensure the precision presented in Tables 1, 2. 

III. COMPARISON WITH EXACT-DIAGONALIZATION DATA AND DISCUSSION 

For the pure J\ = J2 Heisenberg model, the series for the ground-state energy E and on-site magnetization to read 
E/2N = -0.33521 + 0.00019A 2 + <3(A 3 ) = -0.33502 + 0(A 3 ), (23) 

to = 0.3034 + 0.0022A 2 + 0(A 3 ) = 0.3056 + 0(A 3 ), (24) 

where A = 1. The estimates for E and to are in accord with previous large— S@@ and numerical series results. 
The above series show that in the pure model the second-order 0{V^ M ) corrections are numerically small ( i.e., the 
effective quasiparticle interaction is weak) which means that the theory is self-consistent up to this order. As a matter 
of fact, the corrections AE^ do not exceed 3.5% of the MFT values in the whole range —1.2 < a -w30 where the 
on-site magnetization to is positive. Notice that the second-order correction AE^> becomes negativetj for a < —0.8 
and a > 1.4, Table 1. The second-order correction to the on-site magnetization Ato/ 2 ) is less than 10% as compared 
to too in the interval —0.9 < a < 3.4. AmA 2 ' is negative for a < —0.8 and a > 1.6, Table 2. 

In Fig. 3 we show the SWT results for the ground-state energy of Eq.(|l|) compared with ED data . The ED data 
are slightly size dependent, so that, for sake of clarity, we present only the results for a cluster with TV = 26 spins. 
In the parameter range —0.4 < a < 2.0 we find excellent agreement with ED data. Remember that a s w 2.56 
is the numerical-series estimate for the order-disorder critical pointtl, and a = —1/3 is the LSWT estimate for 
the second critical point of Eq.(|l|). For a > 2.0 and a < —0.4 one observes an increasing discrepancy with ED 
data. At the same time, the theory implies very small corrections to the zeroth-order energy Eq in the whole cange 
— 1.2 < a < 4.0. Analogous behavior of the SWT scries has recently been indicated by Chubukov and Morilld in 
the two-layer Heisenberg antiferromagnet when approaching the quantum paramagnetic phase. It was argued, using 
an appropriate bosonization schemea, that the spin-wave description fails due to the enhanced longitudinal spin 
fluctuations near the quantum paramagnetic phase. Near the classical transition point a c = —1/3, where the melting 
of the Neel state is a result both of the frustration and zero-point spin fluctuations, the discussed discrepancy is less 
pronounced. 

As a whole, the second-order correction to the on-site magnetization ArrS 2 \ Fig. 4, improves the estimates for the 
helicoidal a c ph —1.2 and paramagnetic a s ~ 5.0 critical points. However, the estimate for a s is still twice larger than 
the numerical-series result a s « 2.56. At the same time, for instance at a = 2.5, the correction Ato' 2 ' is less than 3% 
from the mean-field value too. In Fig. 4 we also present extrapolated ED data for the on-site magnetization deduced 
from the antiferromagnetic structure factor M 2 = ^ ,.(— l) l+J < SjSj > using the cluster set N — 16,18,20,26 
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and the extrapolation ansatz M| = (oo) + constN~^ . JChe ED estimate m ~ 0.28 at a = 1 is smaller (but in 
agreement with other ED data using also the iV = 36 clustered) than the SWT and numerical-series results mentioned 
above. The above discrepancy probably indicates the importance of the next-order term in the extrapolation ansatzeS. 
Notice that our ED estimate a s k, 2.4 is in reasonable agreement with the numerical-series result a s « 2.56. 

The SWT result a c « —1.2, Fig. 4, also exceeds the estimate obtained through an exact diagonalization of small 
clusters with N — 16, 18, 20 and 26 spins. The ED data predict a phase-transition point a c which is close to the 
classical result a c = —1/3. Here the situation is analogous to the onc.iu.Jiic J\ — J2 model containing frustrating 
diagonal bonds. For this model, both LSWT and the recent ED dataQtj'E3 predict a complete melting of the Neel 
phase at J2/J1 ~ 0.4 , whereas an inclusion of the spin-wave interaction changes the estimate at least up to the 
classical transition. point J2/J1 = 0.5ESE3. There is a smaller SWT estimate J2/J1 ~ 0.43 based on the standard 
large-S expansion^, however it comes from large negative corrections Am' 2 ' reducing the effect of the diverging Am' 1 ' 
term, Eq.(||). 

In conclusion, we have studied the influence of quantum spin fluctuations on the zero-temperature phase diagram 
of a S = 1/2 square- lattice Heisenberg antiferromagnet with two kinds of nearest- neighbor exchange bonds. The 
phase diagram is found to contain a helicoidal magnetic phase and a finite-gap quantum paramagnetic phase, both of 
them being separated through continuous transitions from the classical Neel phase. We constructed an appropriate 
spin-wave expansion and used exact-diagonalization data to find the changes of the phase diagram implied by the 
zero-point spin fluctuations. The latter approach reveals in particular the limits of the spin-wave description coming 
from the enhanced longitudinal spin fluctuations near the phase boundaries. The disagreement with. ED data is 
strongly revealed in approaching the quantum-paramagnetic-phase boundary. In a forthcoming papered we present, 
in particular, ED-based analysis concerning the finite-gap quantum paramagnetic phase. 
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APPENDIX A: DYSON-MALEEV VERTICES 

The normal-ordered Dyson-Maleev interaction Vdm can be written in the following for 

^=-|fE a ( i + 2 - 3 - 4 ) 

1-4 

X 0.30.4 

+ vWpl&c&Pi + V {7) o\o\f3lf3\ + F (8 )/3i/3 2 a 3 a 4 + y (9) /?I/3^ 2 /?i] , (Al) 

where the sum runs over the vectors ki,k2,k3,k4 in the magnetic Brillouin zone and the Kronecker function A(k) 
expresses the conservation of momentum to within a reciprocal-lattice vector G. For the Hamiltonian the nine 
complex vertex functions are 

V«(1234) = W(1234)(a:i [a£(7i-4 - £371-3-4) - (71 - £371-3)] 

+ £ 2 [£4(72-4 - £372-3-4) - (72 ~ £372-3)] ) , (A2) 

V^(123A) = 2U(12M)(x 1 x* 2 [^(^ 7l _3-4 - 71-4) - (^7i-3 -71)] 

+ [x\{xl 72-3-4 - 72-4) - (x* 3 72-3 - 72)] ) , (A3) 

V* 3 >(1234) = 2W(1234)(xi [^(71-3-4 - £371-4) - (71-3 - £371)] 

+ £2 [£4(72-3-4 - £372-4) - (72-3 ~ £372)] J , (A4) 

V (4) (1234) = 4U{123A)(x lX * 2 [14(^71-3-71) - (£371-3-4-71-4)] 

+ [£4(£372-4 - 72) - (£^72-3-4 - 72-4)] ) , (A5) 

V^(123A) = 2U(12M)(x* 2 [(x* 7l _ 3 -4 - 71-4) - £4(^71-3 - 7i)] 

+ x l [(£3 72-3-4 - 72-4) - £4(£372-3 - 72)] ) , (A6) 

F^(1234) = 2W(1234)(^x 2 [z 4 (£ 3 72 - 72-3) - (£372-4 - 72-3-4)] 
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+ [£4(£37i - 71-3) - (2:371-4 - 71-3-4)] J , 
W 7 )(1234) = W(1234)(a; 1 [0=4(71-3 - z 3 7i) " (71-3-4 - ^71-4)] 

+ £2 [£4(72-3 - £372) - (72-3-4 - £372-4)] J , 
VW(1234) = W(1234) (x* 2 [(£371-3 - 7l ) - ^(x* 7l _ 3 _ 4 - 71-4)] 

+ X\ [(^372-3 - 72) - 4(^372-3-4 - 72-4)] J , 

^(1234) = W(1234)(a:3 M71-3 - £371) - (71-3-4 - £371-4)] 
+ x\ [^4(72-3 - £372) - (72-3-4 - £372-4)] J , 

-Uk/wk are declined by Eq.flll 



where ^(1234) = U1U2U3U4 and the Bogoliubov coefficients and x^ 
complex lattice structure factor given by 

7k = - no cos k x + ino sin k x , 



(A7) 

(A8) 

(A9) 

(A10) 
7k is a 

(All) 



where = (cosk x + cosk y )/2, Kq = (1 — a)/4. | _In the symmetric Ji = J 2 model, the vertex functions are real and 
reduce to the well-known Dyson-Maleev verticesEj. 

APPENDIX B: CALCULATION OF AM (2) 

A straightforward way to find the second-order correction to on-site magnetization Am' 2 ' is to introduce in H a 
perturbing staggered magnetic field h (see, e.g., Ref. |l6| ). The perturbed Hamiltonian reads 



H(h) =H- hNfn = E (h) + H Q (h) + V DM , 
where the on-site magnetization operator has the form 



= s--y 



l u k| 



1 K \ \ 

(4«k + PIP*.) - 2ul{x k alf3l + h.c.) 

k £k 



N 



In the new basis of Hq(H) the second-order correction to the ground-state energy reads 

AE {2] (h) = < ° h \ VDM \ ph ><C P h \ Vr > M \° h > 

[ph] 

where \ph > are four-particle eigenstates of Ho(h). 



(Bl) 
(B2) 

(B3) 



Am' 2 ' can be obtained by applying the Hellmann-Feynman theorem to Eq.(Bl) and then using Eq.(B3) 



N oh 



(B4) 



There are two contrib utio ns to Eq.(p3) Ami 2 ' and Ato[ 2 \ coming respectively from the 0(h) terms in the dominator 



and numerator of Eq.(B3). A simple calculation using 



E k (h) =£ k + (1-k)— + 0(h 2 ), 

£k 



(B5) 



gives the correction Ami 2 , Eq. (I 



The second contribution Am[ ' results from the O(h) corrections in the numerator of Eq.(B3). It is easy to check 



that the field corrections to the four-particle states \p > do not contribute up to the required 0(h) order, so that it is 
enough to put in Eq.(jB3|) the unperturbed state \p >, and to use the following expression for the perturbed ground 
state \0h > 



\0 h> =\0>-2hJ2^§r\k>+O(h 2 ) 



(B6) 



Here |k >= o^/^O >. Then a straightforward calculation of the respective matrix elements yields Eq.(|21|) for A 



,(2) 
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FIG. 1. The classical helicoidal state in a square-lattice Heisenberg antiferromagnet with two kinds of regularly distributed 
nearest-neighbor exchange bonds, Ji and J2. The helicoidal phase is stable for a = J2/J1 < —1/3. The spin orientations at A 
and B lattice sites are defined by the angles 9„ — n$ and 9 n = n$ + n, respectively, where n = 0, 1, 2, .... The state is shown 
for $ = tt/12, i.e., a = -1/(1 + y/3) , and n = 0, 1, .., 7. 

FIG. 2. Phase diagram of the S — 1/2 system in a linear spin- wave approximation (LSWT). The solid lines represent the 
on-site magnetization m as a function of the parameter a = J2/ Ji . The dashed line shows the diverging first-order correction . 
There is a small region around the classical transition point a c = —1/3 (—0.336 < a < —0.331) where the magnetic states are 
melt. For a > a s a finite-gap quantum paramagnetic phase is stable. a s w 10.6 in LSWT. 

FIG. 3. Ground-state energy of the model vs a = J2/J1, S = 1/2. The open circles are ED data for N = 26. ED data for 
N — 16, 18, 20 are not shown but the deviations from the N — 26 cluster are very small. Eo is the zeroth-order mean-field 
result. 

FIG. 4. On-site magnetization in the Neel phase vs a = J2/J1 calculated up to second order in Vdm, S = 1/2. mo is the 
zeroth-order mean-field result and Am' 2 ' is the second-order 0(V^ M ) correction. The stars represent extrapolated ED data 
for the on-site magnetization deduced from the antiferromagnetic structure factor Ml = ^\ (— l) l+J < SiSj > using the 

cluster set N = 16, 18, 20, 26 and the extrapolation ansatz M s 2 = M 2 (oo) + constN~i 



TABLE I. Data for the ground-state energy E of the Neel phase calculated up to second-order corrections AE^ in respect 
to Vdm- a — J2/J1, E = Eo + AE^K Eo is the zeroth-order mean-field energy. S = 1/2. 



a 


E /2N 


A£ (2) /2AT 


E/2N 


a 


E /2N 


AE (2) /2N 


E/2N 


-2.0 


-0.21563 


-0.036 


-0.252 


0.4 


-0.29085 


0.00005 


-0.29080 


-1.8 


-0.21811 


-0.023 


-0.241 


0.6 


-0.30429 


0.00009 


-0.30420 


-1.6 


-0.22094 


-0.014 


-0.235 


0.8 


-0.31909 


0.00015 


-0.31894 


-1.4 


-0.22421 


-0.0078 


-0.2320 


1.0 


-0.33521 


0.00019 


-0.33502 


-1.3 


-0.22604 


-0.0054 


-0.2315 


1.2 


-0.35259 


0.00016 


-0.35243 


-1.2 


-0.22800 


-0.0035 


-0.2315 


1.4 


-0.37117 


0.00006 


-0.37111 


-1.1 


-0.23012 


-0.0021 


-0.2322 


1.6 


-0.39089 


-0.00015 


-0.39104 


-1.0 


-0.23241 


-0.0010 


-0.2334 


1.8 


-0.41168 


-0.00046 


-0.41214 


-0.8 


-0.23755 


0.0002 


-0.2374 


2.0 


-0.43348 


-0.00088 


-0.43436 


-0.6 


-0.24356 


0.0006 


-0.2430 


2.4 


-0.47987 


-0.00206 


-0.48193 


-0.4 


-0.25056 


0.0006 


-0.2500 


2.8 


-0.52964 


-0.00372 


-0.53336 


-0.2 


-0.25870 


0.00037 


-0.25833 


3.6 


-0.63800 


-0.00852 


-0.64652 


0.0 


-0.26808 


0.00017 


-0.26791 


4.0 


-0.69606 


-0.01173 


-0.70779 


0.2 


-0.27878 


0.00006 


-0.27872 
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TABLE II. Data for the on-site magnetization m of the Neel phase calculated up to second-order corrections Am' 2 ' in 
respect to Vdm- a = J2/J1, m = mo + Am' 2 ', mo is the zeroth-order on-site magnetization. S = 1/2. 



a 


m 


Am' 2 ' 


m 


a 


m 


Am' 2 ' 


m 


-1.2 


0.1726 


-0.14 


0.03 


1.4 


0.2996 


0.0016 


0.3012 


-1.1 


0.1819 


-0.08 


0.10 


1.6 


0.2952 


0.0009 


0.2961 


-1.0 


0.1911 


-0.04 


0.15 


1.8 


0.2894 


-0.0002 


0.2892 


-0.8 


0.2095 


0.001 


0.211 


2.0 


0.2824 


-0.0016 


0.2808 


-0.6 


0.2274 


0.012 


0.239 


2.2 


0.2743 


-0.0034 


0.2709 


-0.4 


0.2444 


0.0112 


0.2556 


2.4 


0.2652 


-0.0054 


0.2598 


-0.2 


0.2599 


0.0075 


0.2674 


2.6 


0.2553 


-0.008 


0.247 


0.0 


0.2734 


0.0044 


0.2778 


2.8 


0.2446 


-0.010 


0.235 


0.2 


0.2846 


0.0026 


0.2872 


3.0 


0.2331 


-0.014 


0.219 


0.4 


0.2931 


0.0019 


0.2950 


3.2 


0.2210 


-0.017 


0.204 


0.6 


0.2990 


0.0019 


0.3009 


3.4 


0.2081 


-0.021 


0.187 


0.8 


0.3024 


0.0021 


0.3045 


3.6 


0.1946 


-0.025 


0.170 


1.0 


0.3034 


0.0022 


0.3056 


3.8 


0.1805 


-0.030 


0.150 


1.2 


0.3024 


0.0020 


0.3044 


4.0 


0.1657 


-0.036 


0.130 
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